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Abstract 


In this work, we investigate inflationary cosmology using scalar field theory deformed by the 
generalized uncertainty principle (GUP) containing a linear momentum term. Apart from being 
consistent with the existence of a minimum measurable length scale, this GUP is also consistent 
with doubly special relativity and hence with the existence of maximum measurable momentum. 
We use this deformed scalar field theory to analyze the tensor and scalar mode equations in a de 
Sitter background, and to calculate modifications to the tensor-to-scalar ratio. Finally, we compare 
our results for the tensor-to-scalar ratio with the Planck data to constrain the minimum length 
parameter in the GUP. 
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I. INTRODUCTION 


Even though we do not have a full theory of quantum gravity, it is expected that various 
approaches to quantum gravity may be used to gain some phenomenological information 
about quantum gravitational processes. One common feature among most of these ap¬ 
proaches, including string theory and semi-classical black hole physics, is the existence of a 
minimal observable length [l-7|. The existence of a minimal length leads to the modihcation 
of the Heisenberg uncertainty principle to a Generalized Uncertainty Principle (GUP). Since 
the Heisenberg algebra is closely related to the Heisenberg uncertainty principle, the GUP 
deforms the position-momentum commutation relation by an additional quadratic term in 


momentum 


u, 


[x,p] = i(l-f/3p^), (1) 

where parametrizes the minimum length. Phenomenological aspects of GUP effects 
have been analyzed in several contexts such as the self-complete character of gravity j^. 


the conjectured black hole productions at terascale 
oscillations [ll|. 


in| and the modihcations of neutrino 


Recently, another interesting form of the GUP has been proposed in 12 h 14| that contains 
a linear term in momentum 


[xi,pj] = i 




a p5ij + 


PiPj 

P 


{p^6ij 3piPj) 


( 2 ) 


where a parametrizes the minimum length. This form of GUP is consistent with doubly 
special relativity (DSR) and black hole physics. DSR [^, 161 is an extension to special 
relativity in which there is an observer independent maximum energy scale, usually the 
Planck energy, apart from the observer independent maximum velocity i.e., the velocity of 
light. Since this GUP is consistent with DSR, this GUP implies the existence of a maximum 
observable momentum, apart from the existence of a minimum measurable length 12l-ll4|. 

Based on some dimensional arguments, one would expect that the fundamental length 
scale would be the Planck scale, but this is not necessarily true. In string theory for example, 
the st ring length scale could be one or two orders of magnitude greater than the Planck length 
scale 17j, or even more in the ADD model 18|, ll9|]. This is why the minimum length scale 
from the GUP is mrametrized by the free parameter a. Upper bounds on a have been 
calculated in [3, 20|, and it was proposed that the GUP may introduce an intermediate 
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length scale between Planck scale and electro-weak scale. Recent proposals suggested that 


these bounds could be measured using quantum optics techniques 
techniques 


22 


2l| or gravitational wave 


23| . If this happens, then it can be considered a milestone in quantum 


gravity phenomenology. In a series of papers, various phenomenological implications of the 
new model of GUP were investigated 


24l-l31|. A detailed review along the mentionec 


of minimal length theories and quantum gravity phenomenology can be found in 


32 


lines 


Ml- 


One possible area in which minimum length could yield observable signatures is inflation. 
In inflation, quantum fluctuations in the early universe gets stretched out to cosmolo gica l 


sizes, thus providing an explanations for structure formation from first principles 


35 


Inflationary models require a minimum of 60 e-folds of inflationary expansion to explain 
cosmological puzzles. This means that inflation expanded Planckian scales to astrophys- 
ical scales, and hence could act as a microscope for probing short distance scales. Thus, 
suggesting that Planck scales could have observable effects on the universe today. 

However, there is the possibility that Planck scale effects gets diluted with the expansion 
of the universe. This is the case in black hole physics, where Hawking radiation does not 
carry imprints of trans-Planckian physics 39|, l40|]. In cosmology however, Brandenberger 


and Martin 


41| considered various models of trans-Planckian physics and found that some 


of those models yield some late time consequences while others do not. These models were 
based on modified dispersion relations. 

Another approach is the model suggested in [42, l43[ based on GUP ([I]), and was used to 
investigate the effect of minimum length on inflation. Since the GUP is independent of an 
exact theory for quantum gravity, this approach remains valid for many different theories 
for quantum gravity. The numerical solution for tensor perturbations in this model has 


been studied in 


44 


45| . and it was observed that the effect on the power spectrum is of the 


order ^/J3H. It has also been demonstrated in ^ that the presence of cutoff affects the 
boundary terms in addition to the bulk terms in the Lagrangian. However, all this analysis 
was done using the GUP in Eq.([T]). In this paper, we use the linear GUP ([2]) for performing 
a similar analysis. We calculate both the scalar and tensor perturbations, in addition to the 
tensor-to-scalar ratio for a field theory deformed by the linear GUP in de Sitter background. 
We compare our results for the tensor-to-scalar ratio with the Planck data and use it to rule 
out certain ranges for the parameter a. 
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II. SCALAR AND TENSOR FLUCTUATIONS WITH MINIMAL LENGTH 


In this section, we investigate how the action of a scalar field theory is deformed by the 


GUP. The action of scalar perturbations can be written as 

3 

di) 


46 


T7| 


= 9 / d7](fyz‘^ , 


( 3 ) 


i=l 


where 0 is the gauge-invariant intrinsic curvature, 2 ; = V'o homogeneous 

part of the inflation field. An alternative way to write the action is 


= 


drjd'^y ( {d^uY — S^^diudju -\ - u" 


( 4 ) 


where u = —z(j). Those two actions are equivalent in absence of minimal length up to a 
boundary term. 


Similarly, the action for the tensor perturbations could take the form 46l. l47 1 

d7jd?ya?{rj)d^K‘ jd^hd ■, (5) 


.(1) _ Ml 

~ QAtt 


or 


42 ) = T / dyd^y ( dyP^d^P^j - d^^drPi^dsPY + -PdP\ 


J i 5 


where hij is the transverse traceless part of tensor perturbations of the metric, and 


= 


— Pa{ri)h](y). 
327r 


( 6 ) 


( 7 ) 


In the presence of minimal length, the equations of motion were derived in 


42 


46| using 


the GUP ([T]). Here we review the derivation of the equation of motion from 4^^ but using 
the linear GUP in Eq.(l2]). 

Since the GUP implies minimal length, we transform to proper space coordinates x* 
instead of the co-moving coordinates ?/* 


S = dr] d^x 


2a3 




~ 


3a' 


2=1 




( 8 ) 


2 = 1 


where the prime denotes derivative with respect to rj. Now, by dehning 

(01,02) := j d^x 0t(x)02(x) 

x*0(x) := x''(j){x) 
p^(f){x) := -id^i(t){x), 
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(9) 

( 10 ) 

( 11 ) 
















we can write the action in the form 


s = j ^Hv)Mv)4>) - (0,pV) } , 

where the operator A{ri) is dehned as: 

A^(a„ + i- Vpv-it”' 

V “ S' ‘ 


( 12 ) 


(13) 


In Eq. ffT^ . the helds are time dependent abstract vectors in a Hilbert space representation 
of the commutation relations 

[x\p^]=iS^^. (14) 

We aim here to incorporate the GUP by modifying the underlying three-dimensional 
position-momentum commutation relations in Eq.(IT4j) with that in Eq. ([2]). A convenient 
Hilbert space representation of the new commutation relations is on helds 0(p) over auxiliary 
variables p® 

x^(j){p) = idpi(j){p), (15) 


p®0(p) = 


P 


with scalar product 


1 -|- ap — 




(01,02)=/ d p 0t(p)02(p). 


(16) 


(17) 


The action, as given in its abstract form, Eqs.([8]) and (IE]), with the new commutation 
relation in Eq.([2|) underlying, can now be written in the p-representation as 


'“’'I 


ip^<p- 


f) _ ^ 

^ a 1 -I- ap — 2a^p^ ^ a 


aVl0l^ 


(1 -I- ap — 2a2p^)^ 

(18) 

The p modes are coupled, so we dehne a new coordinate system as follows 


r) = p, /c® = ap® exp(ap — a"^p ), 


in which the k modes decouple, because 


d--d - ^ P 

" al + ap-2aV^"‘ 




(19) 


( 20 ) 


With a few calculations, we obtain 


p^dp = k^dk 


^-S(ap-a^ p^) 


a? {1 + ap — 2a^p^) ’ 


( 21 ) 
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and the action now reads 
^ = [ dri [ 


' <a^ /ef3 


(rk —-u 
2a^ 


^r, - 3- ) 0fc(r/) 


hl0fc(h)P 


with u and /i dehned by 


u = 


g-3(op-o2p2) 

1 + ap — 2 a^p^ 


h = 


aV' 


(1 + ap — 2a‘^p‘^y 


( 22 ) 


(23) 


This would yield the modihed equation of motion with new functions u and p as follows 


46| 




z 

z 


z' u' 


Z V 


= 0> 


(24) 


where we dehned v^. = a^4>j^. 


It may be noted that the form of the functions p and u is different from their form for 


the regular GUP l4^. Since these functions are important to determine the effect on 


GUP on inhation, we expect that the effect of the linear GUP on inhaton will be different 
from that of the quadratic GUP. However, apart from the form of p and the remaining 
calculations are very similar. So, we could just use the expressions used for regular GUP 


46| 


'^k '^k + ( /^ 


z' iz' 


V 


Z V 
J! 




< + -<+ h--)^^ = 0 


// ^ / / 

Ur H Ur + p 


V 


a" 

a 


a! v' 


// ^ 

“fc + j, “-fc 


-Ur 


a V 


h- 


“fc = 0 


= 0 


(derived from 

(derived from 

(derived from 

(derived from 
:( 2 ) 


(25) 

(26) 

(27) 

(28) 


In this paper, we use Eq. fl26|) for scalar perturbations, since Sg is more commonly 
used in the literature because of its similarity with a massive free scalar held in Minkowski 
space-time j^. For the tensor perturbations, we use Eq. (HT} as was used in 44, 47|. We 
consider the case of de Sitter space in which a = —1/Hp, z' jz = a'/a, and z” jz = a!' j a. 


III. TENSOR PERTURBATIONS 


In this section, we analyze tensor perturbations deformed by linear GUP. The equation 
of motion for tensor perturbations can be written as Eq. (ETj), 


V 




u 


a 

a 


a' v' 


a V 


“fc = 0 


(29) 
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We begin by expressing v and /i in terms of k which was defined in Eq. ffT^ . To do that, we 


introduce the function V, which we define as the inverse of xe ^ ^ such that 

V(xe-^-^') = X. 

From Eq. ([19]), this allows us to write p in terms of k as 

p = 

a 

where we dehned w = —ak/a. Thus, z/ and p could be written as 


^3[V(w)+V^iw)) 


V = 


P = 


a^V^{w) 


and 


1-V{w) -2V^{w)' "" a^l-V{w)-2V^{w)) 

v' a’ V{w) (4 + 7V{w) - 12V\w) - 12V^{w)) 


2 ’ 


(30) 


(31) 


(32) 


(33) 


V a 1 - V(w) - 2V^w) 

The equation of motion has a singularity at tc = —1, because V{—1) = —1 and the 

denominator 1 — l/(—1) — 2V‘^{—1) = 0. This corresponds to conformal time rj = where 
Pk is when tensor perturbations with different co-moving momentum k reach the cutoff 
p = 1/a. For de Sitter space p = —1/Ha, and hence pk = —1/akH. 

The V function has a series expansion near w = —1 (see appendix) 


T./ N 3 1 2 209 3 

V(w)^-l+p--p + —p 


(34) 


with p — + 1). By repeating the analysis done for the regular GUP in we extract 

the most singular terms of the equation of motion to get the leading behavior of 




1 

2y 






(35) 


y \6cr^ 

where p = pk{l — y), a = aH, and the dot denotes derivative with respect to y. The 
parameter a is the ratio between the minimum length scale a and the Hubble length scale 
H~^, and will be useful because in subsequent equations, a and H are always multiplied by 
each other. 

It may be noted that ignoring the term, we obtain 


u'/ + uluk = 0, 


where ul = A/ply, and 


A = 


6 ( 7 ^ 


1 

2 ' 


(36) 


(37) 
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The solution is approximated by the WKB form 


1 n 

J UJk{v')dv') 

This choice of vacuum is called the Bunch-Davies vacuum, and we get 

^ (4^1) 


( 38 ) 


(39) 


Even though we cannot ignore the term, it is possible to use the solution obtained 
in Eq. (15^ to obtain a deformed solution for the actual equation. Thus, we modify the 
prefactor of the exponential to obtain an exact solution for Eq.( 


F{y) = ^exp{-2i^/Ay). 


The general solution for Eq. fl3^ can now be expressed as 


(40) 


Ui(y) = C\F(y) + C.F-(y), 


where the constants are constrained by the Wronskian condition 42l. |44 1 


(41) 




(42) 


leading to 


i2 1^12 _ _ 2\/M / 1 


\c+Y-\c.Y = 


A. «= —1^+’' 


(43) 


It may be noted that this is y dependent, but this is acceptable since we only need an 
approximate solution for Eq. fl29p near the singularity to solve the equation numerically. 
Comparison with the Bunch-Davies solution suggests that C_ = 0, which was confirmed in 


4i,l47|. 


The solution fl4(I]) is an approximate solution to Eq. fl2^ : we can get a more accurate 
solution by extracting the subleading behavior of Eq. fl29p using the method of dominant 
balance. Thus, using uj,{y) = F(|/)(l -|- ei(|/)), we obtain the equation 


which has the solution 


1 . .20 /3kll 
2y 27 V 2 2 / 


/ 3 20 2A 


(44) 


( 45 ) 













abs(M/;) 



FIG. 1: This figure shows how each mode evolves as a function of rj for a = 0.01. IC-I was taken 
to equal zero. 


Now we can mprove the solution by u{y) = F{y){l + ei( 2 /))(l + € 2 ( 1 /)), to obtain 
... 1 _ 1 f 40\/6^ _ 8 /3 , , 400 /I 8 . . 8 /3 


+ 27 


27V2'^ + i^VT'^+9V2 '■ 


(46) 


The solution to this equation can be written as 

72\/6 


(^2{y) = 


243 


(1 + 3A)y^/^. 


(47) 


Thus, the analytic solution near the singular point is given by 

“i = C+F{y){l + Ci{y))(l + £ 2 ( 9 )) + C'_F'(!/)(1 + £i(!/))(l + £ 2 ( 9 )) 


(48) 


This solves the equation of motion near the singular point when the mode is well inside 
the horizon. To find the behavior of the tensor perturbations when the mode is well outside 
the horizon, we solve the equation of motion numerically. This is done by evolving the 
analytic solution (I48|) when rj ^ rjk numerically. In Fig{Tl we plotted the absolute value 
of for a = 0 . 01 , and we see that it starts at a specihc value for 77 , then oscillates and 
increases until it goes to inhnity as rj goes to zero at present time. 

The range of possible values for a = aH is a < 1. This is because the lowest bound found 
on the GUP parameter a is a < 10® from the anomalous magnetic moment of the muon 


48| . In addition, during inflation, the universe expands exponentially such that a cx 


and most inflationary models require a minimum of 60 e-folds of inflationary expansion, i.e. 
Ht > 60, and that inflation ends after about 10 ®‘^s, or 10^° in Planck units. This means 
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FIG. 2: The dependence of the tensor power 
independent of the value of k, and we obtain tl 


,1/2 

T 



1 /2 

spectrum Pj,' on a when C- = 0. The plot is 
1 /2 

le standard result of Prp = 0.159155 when a ^ 0 


that approximately if ~ 10 The bound on a and the value of H during inflation suggest 
that a < 1. Also, a = 1 means the minimum length scale a equals the Hubble length scale 
so its natural to investigate a < 1. 

The condition for horizon crossing in terms of k can be found from the relation between 
k and p 

F = (49) 


At the horizon, we can express p = H. Thus, in de Sitter space, we obtain 


^horizon 


1 

1 


(50) 


We are interested in the asymptotic values of |n^|, when p/H —)■ 0, which we implemented 
in the numerical solution as p/if —)■ 0.01 


Pasymp 


0.01 


exp(a/100-aVl0000). 


The tensor power spectrum in near de Sitter space can be expressed as 


44| 


pl/2 

r^T 


/ 

U'k 

/ 27r^ 

a 


(51) 


(52) 


This is evaluated at the asymptotic value Pasymp- In FigJ21 we plotted the dependence of the 
tensor power spectrum on a, when C- = 0. The plot is independent of k, which was checked 
by using different values for k and getting the same results. We obtain the standard result 
of = 0.159155 when a ^ 0. 
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IV. SCALAR PERTURBATIONS 


In this section, we analyze scalar perturbations deformed by linear GUP. The equation 
of motion for scalar perturbations fl2B]) can be written as 


y . / a" 

+ I A*- — I *>1 = 0. 


(53) 


where in de Sitter space z”/z = a” ja. 

Following the same steps as for tensor perturbations, we extract the most singular terms 
as 




which has the exact solution 


where we dehned 


Vb 

2 


G{y) = 


B = 


6 ( 7 ^ 


Thus, the general solution around the singularity is 


(55) 


( 66 ) 


Vi = DMv) + D.G’(y), 


(57) 


with the Wronskian constraint 


|T>+|2 - |L)_p = 




-VkV- 


(58) 


To hnd an approximate solution to Eq. flS^ . we use the method of dominant balance by 
substituting = G(|/)(l + ^i(|/))(l + ^ 2 ( 1 /)) in Eq. (l5^ to get the two equations 


and 


y 1 5 .20 3 r-\ 

H-G = i—\ -VB-, 

^ 27 \ 2 y’ 


A- 1 ^ -400 1 


which have the solutions 


, , , .20 2B 


and 


&(») = {siV6B + lOOiVB'j 


(59) 


(60) 


(61) 


(62) 
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abs(vt) 



FIG. 3: This figure shows how each mode evolves as a function of rj for H = 0.1 and a = 0.1. \D^ 
was taken to equal zero. 




1 /2 

FIG. 4: The dependence of the scalar power spectrum Pg on a when H = 1 and C- = 0. The 
power spectrum is independent of the value of k, and we obtain the standard result of Pg = 
0.159155 when a ^ 0 


Thus, we obtain the following general solution 


= D+G{y){l + ei(|/))(l + 6(l/)) + D.G*{y){l + UvW + 6(1/)). (63) 


We evolve this solution numerically to obtain the behavior of the scalar modes after 
crossing the horizon. The scalar power spectrum is given by 


pl/2 


/ F 


/ 27r2 

z 


(64) 


In Figini we plotted |n^|, and in FiglU we plotted the scalar power spectrum as a function of 
a. Both plots are qualitatively similar to those of tensor perturbations but slightly shifted. 
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V. TENSOR-TO-SCALAR RATIO 


The tensor to scalar ratio is given by 


r = 




(65) 


where At and As are the tensor and scalar amplitudes respectively, and are defined by [49 1 


AT{k) = -P}/^ 


10 


T ’ 


Adk) = 


( 66 ) 


In FigJSl we show the tensor to scalar ratio r as a function of a for \C-\ = \D^ \ = 0. For 
small values of a, r oscillates rapidly about its standard value, and at certain values for a 
it gets much larger than the standard value, but r decreases monotonically for a > 0.2. We 
can see why these rapid oscillations occur by looking at Figs. |2] and H] for and 
respectively. These two figures are similar, but slightly shifted from each other. So at the 
values of a when Pj,' gets close to zero, the ratio r also gets close to zero. But at the values 
of a when Pg gets close to zero, the ratio r increases rapidly. 

The upper bound on r from the Planck 2015 results 50| is r < 0.11. The value r = 0.11 
corresponds to the horizontal line in FigJSl and we see that it excludes some values for a, 
and since a = aH, this constrains the possible values of the minimum length scale a once 
H is determined in a specific inflationary model. 

The tensor to scalar ratio for the GUP ([T]) has been investigated in It was demon¬ 
strated that the value of r oscillates around the standard result, however the size of these 
oscillations tends to increase without a bound. So, the value of r becomes bigger than the 
Planck data bound on it for large values of the GUP parameter. 


VI. CONCLUSIONS 

Inflationary models are quite sensitive to short distance physics, which motivated us to 
use a GUP modified framework for analyzing inflation. The GUP deformation we chose is 
consistent with doubly special relativity. This GUP introduces an additional linear term in 
the Heisenberg algebra, and implies both a minimum length scale and a maximum momen¬ 
tum. We investigated the deformation of a scalar field theory in curved background by this 
GUP. We solved the tensor and scalar mode equations in de Sitter background, calculated 
the modified tensor-to-scalar ratio, and compared our results with the Planck data. 
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FIG. 5: Tensor to scalar ratio r as a function of a for IG-I = \D-\ = 0. 


It may be noted that it is possible that the Higgs held plays a role of the inhaton Ml, l52 |. 
Since the Higgs held is related to the electroweak sector of the standard model of particle 


physics, it is possible that inhation can occur in that sector of the standard model. 


possible to generalize the gauge principle for theories with a minimum length scale 53|, 


t is 


which is done by hrst analyzing the full spacetime deformation caused by the generalized 
uncertainty principle, and then converting all the derivatives in the theory to gauge covariant 
derivatives. Such a generalization for the electroweak sector of standard model has also been 


performed in 


54| . It will be interesting to use this theory to study the Higgs inhation with a 


minimum length scale. It may be noted that motivated by the linear GUP, a non-local gauge 
theory has been constructed Mj- However, even though this theory is non-local, it can be 
ehectively treated as a local theory in the framework of harmonic extension of functions. It 
would be interesting to generalize these results to electroweak sector and then analyze the 
Higgs inhation in this model. 


Appendix: The function V{x) 

We introduced the function V{x), which we dehne as the inverse of , such that 

U(xe-"-"') = X, (A.l) 

or equivalently 

= x. (A.2) 

To hnd series expansion of the function V{x) near any point, one could use the Lagrange 
inversion theorem. However, we need the expansion of the series around x = —1, and V{x) 
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FIG. 6; A plot of V{x). 


is not analytic at this point since 

dV 


V{x) 


dx x{l — V{x) — 2V'^{x)) 


(A.3) 


Following the method used in for the Lambert W function, we put p = \/2{x + l)/3 


in Ve 


-v-v^ _ 


= X, and expand in powers of 1 + 1/ 

|p= - 1 = -1 + |(1 + Vy + 1(1 + 1^)= - 1(1 + V-)* - 1(1 + V')= + ... (A.4) 


Then, this series can be inverted to give 


1 2 I 209 3 


1537 


VM = -l+p-gP +^P 


P^ + .. 


(A.5) 


To evaluate V (x) at any point for the numerical solution, we used Halley’s method. This 
method is used to solve a nonlinear equation f(x) = 0 by the sequence of iterations 

2 /(x„)/'(x„) 




2(f'(Xn))^ - f(Xn)f"(Xn)' 


(A. 6 ) 


This method can be used to hnd successive approximations to x = Ve ^ by 


K+i = K 


He 


-V-y 2 


X 


„_y-y2 


( 1 -H- 2 H 2 ) - 


-a;)(4y3+4y2-5V-2) 


(A.7) 


2(i-y-2y2) 

In Figjni we plot the V function based on this numerical approximation. 
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